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We describe a new code to search for signatures of cosmic strings in cosmic microwave anisotropy 
maps. The code implements the Canny Algorithm, an edge detection algorithm designed to search 
for the lines of large gradients in maps. Such a gradient signature which is coherent in position space 
is produced by cosmic strings via the Kaiser-Stebbins effect. We test the power of our new code to 
set limits on the tension of the cosmic strings by analyzing simulated data with and without cosmic 
strings. We compare maps with a pure Gaussian scale-invariant power spectrum with maps which 
have a contribution of a distribution of cosmic strings obeying a scaling solution. The maps have 
angular scale and angular resolution comparable to what current and future ground-based small- 
scale cosmic microwave anisotropy experiments will achieve. We present tests of the codes, indicate 
the limits on the string tension which could be set with the current code, and describe various 
ways to refine the analysis. Our results indicate that when applied to the data of ongoing cosmic 
microwave experiments such as the South Pole Telescope project, the sensitivity of our method to 
the presence of cosmic strings will be more than an order of magnitude better than the limits from 
existing analyses. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

There has been a recent renaissance of interest in cos- 
mic strings (see e.g jlj ) as a mechanism contributing to 
the spectrum of primordial cosmological perturbations 
and the corresponding anisotropics in cosmic microwave 
anisotropy (CMB) maps. In part this renewed interest is 
sparked by the realization that many supergravity mod- 
els of inflation also lead to the production of strings at the 
end of the period of inflation [3] . In addition, many mod- 
els of inflation in the context of string theory predict the 
formation of a network of strings at the end of the infla- 
tionary phase Q . The conditions under which these cos- 
mic superstrings [3] are stable have been explored in 
Cosmic superstrings may also arise in other approaches 
to superstring cosmology such as the Ekpyrotic scenario 
@ or string gas cosmology 0, [1| ■ 

If matter is described by a particle physics model which 
admits stable strings forming after inflation, then by 
causality it is inevitable that a network of such strings 
will form during the cosmological phase transition which 
corresponds to the symmetry breaking which is responsi- 
ble for the existence of the strings @ . It is inevitable that 
the network of strings contains "infinite" strings (strings 
crossing the entire volume). For applications in cosmol- 
ogy one divides the strings into "loops" (loops of cosmic 
string with a curvature radius smaller than the Hubble 
radius H~^{t), where H{t) is the cosmological expansion 
rate) and "long strings" (infinite strings and loops with 
curvature radius larger than the Hubble radius). The 
causality argument @ (see also p^l - [T2| for some stan- 
dard reviews on cosmic strings and structure formation, 
and [ill for the original references) implies, in fact, that 
at all times after the phase transition a network of long 
strings with a correlation length (mean curvature radius) 
comparable or smaller than the Hubble radius will sur- 



vive. It can be argued [IMlIll 

that a network of (non- 
superconducting) strings will approach a "scaling solu- 
tion" for which at all late times t the correlation length 
of the network of long strings is a fixed fraction of the 
Hubble radius ^. 

Cosmic strings '24'] are described by their mass per 
unit length fi which is usually quoted in terms of the di- 
mensionless number G^, G being Newton's gravitational 
constant. A straight string has a tension which is equal 
in magnitude to /x. Hence, the effective action which 
describes the motion of a straight string is the Nambu- 
Goto action, the same action which describes fundamen- 
tal strings. Since cosmic strings carry energy, they will 
lead to density fluctuations and cosmic microwave back- 
ground (CMB) anisotropics. 

The network of cosmic strings will generate a scale- 
invariant spectrum of cosmological perturbations '25']. 
More relevant to the current paper, cosmic strings gen- 
erate a very specific signature in cosmic microwave back- 
ground anisotropy maps, namely line discontinuities [26| . 
These line discontinuities arise since space perpendicular 
to a cosmic string is a cone with deficit angle given by 

M 

a = SnGfi. (1) 

Since the motion of cosmic strings is relativistic, photons 
passing on different sides of the string moving with a ve- 
locity V perpendicular to the plane spanned by the string 



^ It is likely that the distribution of cosmic string loops also ap- 
proaches a scaling solution; see e.g [lj|4l7| for early simulations 
of cosmic string networks supporting the conclusion that the dis- 
tribution of loops also scales, ilSl for an opposing view based on 
field theory simulations and |1SH2^ for more recent work sup- 
porting scaling for string loops. 
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FIG. 1: Geometry of the Kaiser- Stebbins effect: Photons 
passing on the two sides of the moving cosmic string obtain 
a relative Doppler shift for the observer who is at rest. 



direction and the line of sight between the observer and 
the string will be seen by the observer with a relative 
Doppler shift 
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where j{v) is the relativistic gamma factor (see Figure 
1). Looking in direction of the string, we will see a line 
in the sky across which the CMB temperature jumps by 
the above amount. We will denote this effect, the Kaiser- 
Stebbins effect, as the KS effect in the rest of the paper. 

The conical deformation of space induced by a "long" 
cosmic string has a length of the order the Hubble ra- 
dius, i.e. of the order of t in direction of the string. Since 
cosmic strings are formed in a phase transition and the 
effects of the transition travel with the speed of light, 
the depth of the region affected by the string (in direc- 
tion perpendicular to the string) is t, as shown in [28| . 
Hence, each string between tree and to whose "zone of 
influence" ^ is intersected by the past light cone of the 
observer will induce a line discontinuity in the CMB tem- 
perature map with length gt, where g is a random number 
in the interval < g < 1 which takes into account the 
random angle of the velocity vector of the string relative 
to the plane determined by the string direction and the 
observer's line of sight to the string. The anisotropy pat- 
tern induced by a single string segment contains several 
edges. There is the central line discontinuity ([2]) coming 
from photons which pass on different sides of the string. 
Since the deficit angle in Figure [1] has finite depth (in di- 
rection —v), as explained in the previous paragraph there 
will be two edges (again sharp because the deficit angle 
sharply decreases to zero [l^) with half the value of 6T. 



^ By which we mean the region of space which is deformed due to 
the presence of the string. 



Finally, our modelling of the string network in terms of 
straight segments introduces the sharp "boundary" edges 
perpendicular to the central edge (see also Figure [3]) . 

To detect the line discontinuities in CMB maps pro- 
duced by cosmic strings, it is important to have small an- 
gular resolution. Strings present in the universe between 
the time of last scattering and the present time contribute 
to the signal. However, according to the cosmic string 
scaling solution, the most numerous strings are those 
present at or shortly after last scattering. The Hubble ra- 
dius at that time subtends an angle of about 1.8°. If the 
angular resolution is not much smaller than this angle, 
then the anisotropics produced by these strings will not 
be distinguishable from anisotropics produced by Gaus- 
sian noise with the corresponding coherence length. On 
the other hand, full sky coverage is not essential. Thus, 
ground-based small angular resolution surveys such as 
ACT |2i| or the South Pole Telescope [30] both of which 
have angular resolution of about 1 are ideal to search 
for strings. The Planck satellite experiment [3l| with an 
angular resolution of about 5 will also yield a good data 
set to use, in particular since the systematic errors in the 
data will likely be smaller. 

The KS effect is a part of the "Integrated Sachs- Wolfe" 
(3^ contribution to CMB anisotropics. The primordial 
cosmological fluctuations produced by strings also con- 
tribute to the regular Sachs- Wolfe effect. In contrast to 
cosmological fluctuations produced in inflationary cos- 
mology, those produced by cosmic strings are "incoher- 
ent" and "active" as opposed to "coherent" and "pas- 
sive" [13. The string network is continuously seeding 
growing curvature fluctuations on super-Hubble scales, 
and therefore the fluctuations enter the Hubble radius 
not as standing waves. As a consequence, the angular 
power spectrum of CMB anisotropics does not have [33 - 
ISq the acoustic ringing associated with coherent passive 
fluctuations (37| . Since acoustic ringing has been ob- 
served with recent high precision CMB measurements 
(38l |39| . it is now clear that cosmic strings cannot be 
the main source of cosmological fluctuations. Their con- 
tribution is bounded from the accurate measurements of 
the angular power spectrum of the CMB in the region 
of the first acoustic peak to be less than 10% [i^l which 
corresponds to a value of Gfi of about 3 x 10~^. In the 
literature, one finds slightly stron ger constraints which 
come from pulsar timing data (4ll |42|. However, pulsar 
bounds make use of estimates of the spectrum of gravita- 
tional radiation from cosmic string loops. Since there is 
still a lot of uncertainty about the distribution of cosmic 
string loops, and since the amount of gravitational radia- 
tion from a fixed loop is also rather uncertain, bounds on 
G/z coming from millisecond pulsar timing are not very 
robust. Much more robust signatures come from the long 
and straight strings, signatures which we are discussing 
in this paper. 

Early work to identify the KS signal of cosmic strings 
in CMB anisotropy maps was presented in (43| which 
concluded that the angular resolution of WMAP would 
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not be small enough to resolve the KS signature. After 
the release of the WMAP data 39] , there were two sets 
of analyses introducing new algorithms to look specifi- 
cally for the KS signature. Lo and Wright Q applied a 
matched filtering method, whereas Jeong and Smoot [i^ 
introduced new statistics such as one measuring the con- 
nectedness of neighboring temperature steps or another 
one proposing a decomposition of the temperature map 
into constant, Gaussian and straight string step compo- 
nents. Both groups applied their statistics to the WMAP 
data. From the null results of the searches, a direct upper 
bound on the string tension of G/z < 10~^ could be set, 
a bound not competitive with the existing bounds from 
the matching of the angular power spectrum. 

In [46] , it was proposed to make use of the Canny al- 
gorithm j47(] to search for the KS effect, and the prelimi- 
nary analysis showed that the statistic offers the promise 
to improve the limit on G/i from direct searches for the 
KS signal by a large factor. The Canny algorithm is an 
edge detection algorithm which was previously used in 
image recognition work and metallurgy. It looks for lines 
in a map across which there is a large gradient. Thus, 
the algorithm appears to be well suited to detect the KS 
signature. In this paper we present a new and improved 
implementation of the Canny algorithm and apply it to 
test data. In agreement with [46] (and with the foUowup 
paper [48]), we find the new analysis based on the Canny 
algorithm is able to find or rule out strings with a ten- 
sion greater than an upper bound which is more than 
an order of magnitude smaller than previous bounds de- 
rived from looking for the KS signature of strings. We 
emphasize that the code on which this paper is based 
was developed completely independently from that used 
in [4^ and [1^. It is different in structure and in fact 
is based on a different programming language. The fact 
that the results reported here agree with those in [i^ 
presents a very important check on the methods. 

The Canny algorithm works in position space. Start- 
ing from an anisotropy map, it first produces an "edge 
map" , the edges corresponding to lines in the sky per- 
pendicular to which the gradient of the anisotropy map 
is sufficiently large ^ . The edge detection algorithm must 
be able to take into account the fact that Gaussian noise 
superimposed on top of the cosmic string signal will pro- 
duce large variations in the magnitude of the gradient 
along a cosmic string edge. Given the edge map, a second 
algorithm counts the number of edges of fixed length and 
produces a histogram of edge lengths. Both the edge de- 
tection and the edge counting algorithms contain various 
parameters and thresholds which can be set by the user 
and which can be tuned to give the algorithm maximal 
discriminatory power, and are improved over the original 
code presented in ^]. The values of the thresholds and 



^ As detailed in Section 3, a new aspect of the present code is to 
search for gradients in a range tuned to the expected KS signal. 



parameters will depend on the specific data. The final 
step of the code is a statistical comparison between the 
histograms produced in the previous step for data with 
and without cosmic strings. 

The outline of this paper is as follows: In the following 
section we discuss the construction of the test data maps 
- both the pure Gaussian maps and the maps containing 
cosmic strings. In Section 3 we present the new imple- 
mentation of the Canny algorithm. We first describe the 
algorithm which takes the CMB temperature map and 
converts it into an edge map. Then, we turn to the sep- 
arate algorithm which is used to produce a histogram of 
edge lengths. Section 4 presents the results from the ap- 
plication of the Canny algorithm to the test maps. We 
conclude with a summary of the method and results, fo- 
cusing on refinements of the code which can be made 
to improve the code's discriminatory power, and give a 
preview of future applications. 



II. SIMULATIONS 

In this section we describe the codes which were used 
to create the simulated temperature maps. The maps 
have two components, firstly a Gaussian map with an an- 
gular power spectrum corresponding to the infiationary 
"concordance model" , the second a map of anisotropics 
produced by long cosmic strings according to the KS ef- 
fect. These simulation routines have been created from 
scratch, without using any input from the existing code 
[46j . The simulation routines are written in Interactive 
Data Language (IDL). Both for the inflationary fluctu- 
ations and the cosmic string maps, the theory predicts 
ensemble averages from which particular classical real- 
izations are drawn. To obtain firm predictions, a large 
number of simulations must therefore be run. 

Since we have in mind applying our code to surveys 
with small angular resolution but only partial skycov- 
erage, we work in the "flat sky" approximation [43| in 
which a segment of the sky is approximated by a rect- 
angle. This approximation is made because of computa- 
tional ease, since the basis solutions of the wave equation 
in flat space, the Fourier modes, are much easier to work 
with than the basis functions on a sphere, the spherical 
harmonics. 



A. Gaussian Map 

First we create a square grid temperature map of the 
Gaussian inflationary perturbations. The grid size is set 
by the angular resolution and by the angular size of the 
survey which we want to simulate (the number of grid 
points along an axis is Nmax = L/R, where L is the ex- 
tent of the survey along the axis being considered, and 
R is the angular resolution). The map is constructed 
based on the angular power spectrum of CMB fluctua- 
tions computed using one of the standard codes used in 
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the literature (see below). 

A differential temperature value, AT/T, (correspond- 
ing to the temperature with the monopole subtracted) is 
assigned to each grid point (n, m) (n and m are integers) 
of the map, where each grid point represents a position 
in the window of the observed sky. If L is the length of a 
side of the window of observation in degrees in the first 
coordinate direction, then the angular distance of a grid 
point {n,m) from the edge is x = nL / N^nax- The two 
dimensional vector is designated by x. 

The code is designed to construct a square grid. If 
rectangular areas are needed, we can construct a larger 
square temperature map and cut out the appropriate 
rectangle from this larger grid. The applicable scales are 
less than 60° indicating the applicability of the "flat sky" 
approximation. 

In general, the temperature anisotropy of the CMB is 
expressed in terms of spherical harmonics Yim{9,(j)), 9 
and (j) being latitude and longitude, respectively: 



AT 



(3) 



where aim are the coefhcients of the expansion. However, 
since the observed window will be less than 60 degrees 
(we will take it to be approximately 10 degrees), plane 
waves can be substituted for Yim per the flat sky approx- 
imation: 



AT 



(4) 



(N.B. In the following equations T{x) = AT{x)/T and 
T{k) = AT{k)/T) Comparing these equations, it is obvi- 
ous that the coefficients T{k) correspond to a;,„ so that: 



(T(fc)2) = 



= c, 



l(k)-, 



(5) 



in which the angular brackets stand for ensemble averag- 
ing, and the C; give the angular power spectrum of the 
CMB. According to the ergodic hypothesis, the ensemble 
average is equivalent to the spatial average. Thus, 
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which, as indicated in ([5|) is the square of the width of a 
Gaussian distribution, namely T{k). The probability of 
a given T{k) follows a Gaussian distribution just as the 
aim follow a Gaussian distribution. 

The temperature map in Fourier space, T{k), is ar- 
ranged into a grid in which each grid point is assigned 
physical angular coordinate values kx and ky, which are 
related to the numerical k^n and fcy„ values (each ranging 
between and Nmax — 1) by 



kx {kxn kmax) 

ky — ~^{kyn kmax^ 



(7) 



The value of kmax corresponds to the angular resolution 
of the data set. Note that this setup requires us to have 
an odd number of pixels in our Gaussian map. If an even 
number of pixels is required, the Gaussian map can be 
stripped of a row and column once it is constructed. 
Using 
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in which the wavelength, A is measured in degrees and 
the wavenumber k is measured in inverse degrees, we 
compute the dependence of the degree of the spherical 
harmonic (mode number) l{k) on the magnitude k of k. 
Once this is done, we can look up the corresponding C; 
for each k from the output of a standard CMB simulation 
code (see below). Since l{k) is not necessarily an integer, 
for each value of k, the program then determines the 
value Ci{k) by linear interpolation between the C/ values 
of the two closest values of I corresponding to each l{k). 

Given the Ci values determined as described above, the 
Fourier space temperature map T{k) can be determined 
by 



T{kx,ky) = ^Ci(^k^,k^)/2{gi{kx,ky) + ig2ikx, ky)) (9) 



where gi{kx,ky) and g2{kx,ky) are randomly generated 
numbers obtained from a Gaussian distribution with vari- 
ance one and mean of zero. 

Since the position space temperature map T(x) must 
be real, then T{k) must satisfy T{k) = T*{~k). 

In summary, we derive T{k) from the values of Ci de- 
termined as described above, and T{—k) making use of 
the Hermitian property. Finally, to compute T{x), we 
perform a Fast Fourier Transform on T{k). Figure ^ 
shows the output of a temperature map with pure Gaus- 
sian noise. 

We use the Code for Anisotropics in the Microwave 
Background (CAMB) ^ to generate the CMB tem- 
perature map generated by inflationary perturbations. 
CAMB, like the code on which it is based, CMBFAST 
[Slj . employs the line of sight integration prescription 
to compute CMB anisotropics. The advantage of using 
CAMB is that it computes C;'s for higher Vs. Whereas 
CMBFAST can at best compute the power spectrum to 
I = 3000, corresponding to a resolution of 10 arcmin- 
utes, we use a power spectrum generated by CAMB up 
to / = 22000, corresponding to a resolution of 1.39 ar- 
cminutes. Simulations with higher resolution will enable 
us to place lower bounds on the string tension of cosmic 
strings. Additionally, high resolutions are necessary to 
apply our algorithm in the future high resolution CMB 
experiments such as ACBAR [s^l and the South Pole 
Telescope 

The parameters we use for our simulations are stan- 
dard ACDM concordance parameters based on fitting 
the model to a multi-experiment data set CMBall fssf , 
namely ilth^ = 0.0227, n^h^ = 0.112, = 0.26, 
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FIG. 2: Figure of the gaussian temperature map with a win- 
dow of 10° and a resolution of 1.5 arcminutes. 



Ho = 71.8, Us ~ 0.965, w = —1 (the equation of state 
for dark energy), {l^h^ — 0, Tcmb = 2.726 (the temper- 
ature of the CMB), Yffe = 0.24, and reionization optical 

depth, Treion = 0.093. 



B. Temperature Map from pure Cosmic Strings 

Next we generate a temperature map of superimposed 
straight line segments representing the effects of long cos- 
mic strings obeying the scaling solution. As in the Gaus- 
sian simulation, each grid point is assigned a temperature 
fluctuation AT/T. 

For simplicity we do not include the effects of string 
loops or of the small scale structure on the strings. Ac- 
cording to the current status of cosmic string simulations, 
the effects of long strings dominate the cosmic string 
CMB sky. Small-scale structure on the long strings would 
maintain the KS effect (albeit with slightly reduced am- 
plitude). Thus, in terms of searching for the KS sig- 
nature of strings, the simplifications we use can be jus- 
tified as a first step. Going beyond this approximation 
would require much more involved simulations (like those 
of [l^l) with many more parameters. ^ In future work we 



* Although in principle cosmic string models are determined by a 
single parameter, namely the cosmic string evolution cannot 
be followed either analytically or numerically without making 
additional assumptions, in particular on how cosmic strings in- 
teract. The results obtained in different numerical analyses con- 
cerning the string loop distribution vary by orders of magnitude 
(see e.g. the different results obtained in [TslllTll ). whereas the 
results for the distribution of long strings are comparable. Thus, 



plan to apply our Canny algorithm code to temperature 
maps obtained from more sophisticated simulations that 
include more detailed structure of cosmic strings. 

In contrast to the Gaussian string map which is set up 
in Fourier space and transformed to position space via a 
Fourier transform, the string maps are set up directly in 
position space. As explained in the Introduction, each 
long string whose region of influence intersects the past 
light cone of the window in the sky we are considering 
will contribute a KS signal. We assume that the area 
in the sky affected by the temperature discontinuity of 
a single straight string segment per the KS effect is a 
rectangular box on each side of the string. The depth 
of the box (length perpendicular to the string) is equal 
to the Hubble length [2g|. The size of the box in di- 
rection of the string is taken to be the product of the 
Hubble length and a length coefficient 7. The Hubble 
length corresponds to the time at which the past light 
cone intersects the region of influence of the string. We 
apply the KS signature by adding a positive temperature 
fluctuation to grid points in the box on one side of the 
string, and by adding the corresponding negative tem- 
perature fluctuation to the grid points in the box on the 
other side of the string. 

We are using a toy model for the cosmic string distri- 
bution first introduced in [55] and also used by Pogosian 
et al. in ^40] . It is based on a simple representation of 
the distribution of long strings as described by the scaling 
solution. According to this solution, the distribution of 
long strings is at all times like a random walk with a step 
length which is of the order of the Hubble radius at that 
time. We will denote this length by ^H~^, where 7 is a 
constant factor of order one. Physically, this length rep- 
resents the mean curvature radius of the string network. 
In the toy model, this distribution is approximated by a 
collection of straight string segments of length ^H^^ . 

One can argue that the curvature radius of the network 
of long strings must be of the order of the Hubble radius. 
If the curvature radius is larger than the Hubble radius, 
the strings will sit there as the universe expands until 
the Hubble radius catches up. If the curvature radius is 
less than the Hubble radius, then the strings will oscillate 
wildly, intersect and split off string loops (which then de- 
cay by gravitational radiation) until the curvature radius 
catches up to the Hubble radius. So the curvature radius 
will always be of order of the Hubble radius. The key in- 
gredients in this argument are, firstly, that the effective 
action for a cosmic string is the Nambu-Goto action and 
hence the string dynamics is relativistic, and secondly 
that when strings cross they will intercommute, i.e. ex- 
change ends, thus allowing the production of string loops. 

Since the strings are relativistic, the distribution of 
strings will have changed after one Hubble time. Snap- 



whereas limits based on long strings are robust, those which in- 
clude effects of string loops are not. 
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shots of string distributions at time steps separated by 
Hubble time intervals will, when rescaled to the Hubble 
radius, look like independent realizations of a stochas- 
tic process. Hence, in our toy model of the distribution 
of long strings, we take the string segments to be inde- 
pendently distributed over time intervals greater than a 
Hubble time. 

The input parameters for the string map generating 
function are the number of pixels on each side of the 
window, the number of degrees the window subtends, 
the number of strings per Hubble volume per the scaling 
solution, the string tension, Gfi, and the length coefficient 

7- 

The procedure to simulate cosmic strings is to follow 
the past light cone in time intervals tn-i <t < tn, where 
tn+i — aitn, from the time of last scattering to the 
present time. If ai — tn+i/tn — e then there are ap- 
proximately 15 time intervals. 

Next we project all strings present in one Hubble time 
interval to a fixed time hypersurface at the center of the 
time interval. This is done for all time intervals. On 
the microwave sky, the Hubble length in degrees corre- 
sponding to the comoving distance of the Hubble radius 
in space at the n'th time interval can be shown to be 
(making use of the equation of state and of the Fried- 
mann equation) 

dc{tn+i) = dc{tn)ay\ (10) 

in which dc(t) is in degrees. This recursion relation starts 
with dc{to) = 1.8. 

We need to find all of the strings which for a fixed 
time interval influence the CMB temperature map of the 
observed window. Instead of simulating only strings in 
the observed window, we must consider all strings in an 
extended window. This protects us from missing strings 
that start outside of the observed window but extend into 
it. At each time interval, the extended window is a square 
with each side subtending the degrees of the observed 
window added to twice the Hubble distance for that time 
step. The observed window is the central section with one 
Hubble length on either side. 

Next we compute the number of Hubble volumes the 
extended window subtends at each time step, and the 
number of strings in each extended window. The number 
of Hubble volumes is obtained as the square of the num- 
ber of degrees subtended by each extended window di- 
vided by the degrees subtended by each Hubble distance. 
The number of strings expected in the extended window 
for each time step is the rounded product of the strings 
per Hubble volume (determined as an input parameter 
based on the scaling solution) and the fraction/number 
of Hubble volumes a window subtends for each of the 15 
time intervals. 

Next we loop through each time step and each expected 
string in the extended window to determine where the 
string should be placed in the extended window and to 
place it by constructing a temperature differential. This 
is completed in a separate function in the program. 



To simulate the strings, the program randomly places 
a line segment with a length given by j cosaH^^ , where 
the Hubble length is in radians for the given time 
step. The angle a runs from to 7r/2. Multiplying 
by cos a accounts for a projection from three to two di- 
mensions. Although technically the length of the string 
should be multiplied by a coefficient given by projecting 
a geodesic on a sphere onto a plane, since we are dealing 
with small window sizes and the flat sky approximation 
is appropriate, this complication is unnecessary. The lo- 
cation of the beginning of the string and its direction are 
determined randomly. The location of the beginning of 
the string is given by taking the product of a random 
number between and 1 for both directions on the grid 
and the length of the extended window in radians. The 
direction of the placement of the string, i.e. the angle 9 
of the string, is chosen by a randomly generated number 
between —7t/2 and 7t/2. This is the direction of the string 
as given by the angle from the x-axis. Finally, a binary 
flag is used to determine which side of the string will be 
a positive temperature fluctuation and which side will be 
a negative temperature fluctuation. The value of the flag 
is randomly determined. 

Next the code computes (from the location of the start 
of the string and the direction of the string) the end point 
of the string and each corner of the box around each side 
of the string and each slope for each boundary of the box. 

The routine loops over each pixel in the observed win- 
dow, starting in the lower left hand corner, to determine if 
the pixel lies within the box afi^ected by the cosmic string. 
If it does, the temperature is changed by ST. The sign 
of 6T is determined by the temperature fiag delineated 
above. The temperature fiuctuation is given by 

ST 

— = vriirGfi (11) 

in which T = 2.726 is the background CMB temperature. 
In the above, v represents the root mean square (over all 
strings) value of vj{v), where v is the transverse velocity 
of the string and its relativistic 7 factor is 7(w). Also, 
r a random number between and 1 which adjusts the 
velocity to take into consideration the different velocities 
that the string might have as well as the projection of 
the velocity of the string onto the plane perpendicular to 
the line of sight. Based on recent cosmic string evolution 
simulations we use the value v = 0.15. 

Figure [3] shows the temperature map produced by a 
simulation with a few test cosmic strings (not a scaling 
solution). The temperature boxes produced by the indi- 
vidual strings are clearly visible. Figure U demonstrates 
the corresponding results from a full string simulation 
with A'^ = 1, i.e. one string per Hubble volume. The 
straight line temperature discontinuities produced by in- 
dividual strings are still visible. However, there are a 
lot of overlap regions since a given photon will during 
its trajectory pass close to several strings. The overlap- 
ping problem gets worse for A^ = 10, as clearly visible in 
Figure [S] 





FIG. 3: A cosmic string simulation for a few strings. The light 
areas represent where the temperature fluctuation is positive 
and the dark areas represent where the temperature fluctua- 
tion is negative. 



As a further test of the algorithm described in this 
subsection we show the resulting angular power spectrum 
of temperature anisotropies. The figure (Figure [B]) is for 
a pure cosmic string map with = 10 and G/x = 6 x 
10~®. The range of I values are limited from below by 
the angular size of the simulation box, and from above 
by the angular resolution we have chosen. The error bars 
are standard errors of the mean based on 100 runs. As 
theory predicts [1^, [s^l , the angular power spectrum is 
approximately scale-invariant on large angular scales (the 
deviation from scale-invariance at the lowest values of I 
is presumably a boundary effect), and its amplitude is 
consistent with what is expected. 



C. Sum of Gaussian Map and Cosmic String Map 

Our goal is to test if the Canny algorithm is able to 
pick out the KS signature from strings even if the strings 
are a subdominant component to the CMB fluctuations. 
To test this, we need to produce temperature maps which 
contain both cosmic strings for some value of Gfi and a 
spectrum of Gaussian fluctuations like in the concordance 
ACDM model, except with an amplitude of the Gaussian 
noise which is reduced such that the total angular CMB 
power spectrum remains consistent with either simula- 
tions or observations. 

Here we describe how to find the coefficient, a, for 
each map which will be used to add the pure string map 
to the pure Gaussian map to get a Gaussian map with 



FIG. 4: The CMB temperature anisotropy map produced by 
a scaling cosmic string simulation with A'^ = 1. The disconti- 
nuity lines in the maps produced by the KS effect are clearly 
visible, but there are a lot of overlap regions where a number 
of strings affect the temperature at a fixed point. 




FIG. 5: The CMB temperature anisotropy map produced by 
a scaling cosmic string simulation with TV — 10. In this case 
the effects of overlaps is much more pronounced. 

superimposed strings, 

Ta+sik) = aTG{k)+Ts{k), (12) 

with Tg+s as the temperature of the Gaussian map plus 
pure string map and Tq and T5 as the temperature maps 
of the concordance Gaussian and string simulations re- 



8 




250 



FIG. 6: The angular power spectrum of the CMB anisotropy 
maps of pure cosmic string simulations with values A'' = 10 
and G/i = 6 x 10~*. The horizontal axis is I, the vertical axis 
is l{l + l)Ci. 



spectively. Alternatively, we could use the temperature 
map of the best current data, the five year results of 
WMAP [s^, for Tq. It is important to use consistently 
either the temperature map from experimental data or 
from theoretical simulations for both Tq and Tq^s since 
the Ci's for the two models are not identical. 

We want the model to give the best possible agreement 
with the data, either experimental or simulated. Hence, 
the plan is to adjust the coefficient a so that the combined 
string and Gaussian map, TG+s(fc), fits the data best. 
This means that we adjust the Gaussian map so that 
the corresponding C;'s of Ta+s{k) give an optimal fit to 
the data. Since the error bars on the observed angular 
power spectrum Ci are smallest relative to the signal in 
the range from Imin = 10 to Imax = 220 we shall fit 
the Ci's within this range. Note that this range includes 
most of the first Doppler peak region. Now for values of I 
comparable or larger than that corresponding to the first 
Doppler peak, the contribution from strings is dominated 
not by the Kaiser- Stebbins effect from strings between 
tree and but from string- induced fluctuations at last 
scattering which are not included in our analysis. Thus, 
there is an intrinsic inaccuracy in the determination of 
the value of a. If we were to take Imax to be the value 
where the Kaiser-Stebbins effect ceases to be dominant 
(a value much smaller than 220), then we would obtain 
a smaller value for a and hence a better discriminatory 
power of our algorithm. However, this procedure would 
be worse than the one we have adopted, since we would 
be working with a power spectrum which is a much worse 
fit to the observations in the Doppler peak region than 
the one we are using. Our choice of Imax can thus to be 
considered to be a conservative one. 



The first step is to find the C; values from the pure 
string map. To do this, the routine initially takes the 
inverse FFT of the temperature fluctuation T{x) to com- 
pute T{k). Then, using k = 27rZ/360, the routine com- 
putes the wave number magnitudes k{l) for all values of 
I in the our range. Then the value of Sk corresponding 
to 61 = 1/2 is found. 

Next, the map summation routine finds the wavenum- 
ber magnitude k for all values of k (whose components 
range from —kmax to kmax)- For each value of I in the 
range between Imin and Imax the code finds the compo- 
nents kx and ky of all waves for which the magnitude of 
k lies within 5k of k{l). Knowing fc^„ and kyn we can 
find the corresponding T{k)'^ for each of the wavenum- 
bers. Once all of the T{kY values are found for each 
we find the average to get {T{kY) for each These val- 
ues are the Ci values of the pure cosmic string map for 
all I between Imin and Unax- 

Since 

TG+s{k) ^ aTG{k)+Ts{k) (13) 

we obtain 

{TG+s{kf) - a\TG{kf) + {Ts{kf)+2a{TG{k)Ts{k)). 

(14) 

However, {TG{k)Ts{k)) — because the Gaussian and 
string temperature fluctuations are independent. Hence 

Ci{G+S) = a'^Ci{G) + Ci{S)- (15) 

Now that we have the Ci values for the simulated 
CAMB data, C;(g+S); the Ci's we found for the pure 
cosmic string map, ^(s), and the Ci's for the simulated 
gaussian map from CAMB, Ci(^g)> we can compute for 
each I by viewing (ITSI) as a defining relation for a^. Ob- 
viously, for each value of I we will get a different answer. 
To compute a single coefficient a which best fits all of 
the /'s from Imin to Imax wc perform a linear fit on the 
coefficients, a, for each I. We use the best fit j/— intercept 
of a linear model fit by minimizing the error statistic 
for a. The y-intercept can be used because the slope of 
the best fit to the linear model is negligible. 

Once a is found for a particular string image, the re- 
sult is averaged over all of the generated simulated string 
images to obtain the final value for a which is then used 
for all the string maps of a given Gfi. 

Figure ([7]) shows the resulting CMB anisotropy map in 
a simulation with a large value of Gfi chosen such that the 
strings play an important role. Comparing the map to 
the pure string map of Figure ([5]) we see that the string- 
induced line discontinuities are still visible. However, in 
the case of a lower value of G/i the effects of the strings 
are not visible by eye and we need to resort to a statistical 
analysis to study whether the map is distinguishable from 
that of pure Gaussian noise (Figure dS])). 

We use the algorithms described in this section to gen- 
erate maps with and without cosmic string signals. The 
main question we would like to address now is down to 
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FIG. 7: The smoothed CMB temperature anisotropy map 
produced by a simulation with both Gaussian noise and cos- 
mic strings with Gfi — 3.5 x 10~^ andA'^ = 10. In this case 
the KS discontinuity hues are visible. 




FIG. 8: The CMB temperature anisotropy map produced by a 
simulation with both Gaussian noise and cosmic strings with 
a smaller value of Gfj,, namely Gfi = 2 x 10~^ and A'^ = 10. 
In this case the KS lines are not visible by eye. 



what value of G/i the Canny edge detection algorithm is 
capable of distinguishing between these two sets of maps 
in a statistically significant way. In the following section 
we describe our realization of the Canny algorithm rou- 
tine which turns a microwave temperature map into an 
edge map, the edges standing for lines across which the 



gradient is a local maximum.^ 



III. IMPLEMENTATION OF THE CANNY 
ALGORITHM 

In this section we describe our implementation of the 
Canny algorithm. In its original version [43|, the Canny 
algorithm is intended to find lines in the map with max- 
imal gradients across the line. In this way, the algorithm 
can be applied to the image of a face and returns a map 
where only the pixels of the map with the strongest fea- 
tures are shaded in. Similarly, the algorithm could be 
applied to find crystal defect lines on metallic surfaces. 
Our original work applying the Canny algorithm to CMB 
maps [4g also was based on this idea. 

However, to look for cosmic strings we are less inter- 
ested in the local maxima of the gradient map which over- 
all have maximal amplitude. Rather, we are interested 
in local maxima of the gradient for which the amplitude 
of the gradient is in correspondence to the expected KS 
signal. Thus, in contrast to the original Canny routine 
which uses two thresholds, we here introduce a modified 
algorithm in which three gradient thresholds are made 
use of (as in [11]). The use of the thresholds is described 
below. 

We will first give a brief overview of how the algorithm 
works. First, an optional part of the Canny algorithm is 
to smooth the map in order to eliminate shot noise. The 
second step in the algorithm is to construct a map of the 
temperature gradients. Next, locations must be identi- 
fied which correspond to local maxima of the gradients. 
The maximal average gradient of the pure string map 
determines the three thresholds used. The fourth part 
of the algorithm involves selecting among the grid points 
identified as local maxima those which have the right 
range of magnitudes. Gradients larger than some upper 
threshold are not due to single strings and will hence be 
discarded. Those above a second threshold quite close 
to the expected edge strength will be kept. To take into 
account the fact that Gaussian noise may well decrease 
the amplitude of some pixel along a string edge below 
the above-mentioned threshold, we introduce a third (the 
lowest) threshold and keep the pixels whose gradient is 
larger than that threshold, provided that the pixel is 
connected to a pixel with gradient amplitude above the 
second threshold (and provided that the gradient direc- 
tions are appropriate and that the pixel is connected to 
a pixel in one of the allowed directions relative to the 
gradient). The pixels thus selected form the Canny edge 
map. Given the pixel edge map, a next step in the al- 
gorithm is to identify the edges in the pixel edge map. 



^ In our realization of the algorithm, we impose an additional re- 
quirement, namely that the gradient magnitude is tailored ap- 
propriately to the expected cosmic string signal. 
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The length of each edge is found, and a histogram of edge 
lengths is produced. The analysis part of the algorithm 
then checks if for a fixed value of Gfj, the histograms of 
a pure Gaussian noise map and of a noise map including 
a contribution of cosmic strings (and appropriately re- 
duced amplitude of the Gaussian noise) are statistically 
significantly different. 

In implementing the Canny algorithm, several (three) 
choices must be made beyond fixing the three thresholds 
used. We look for the set of choices which give the best 
differentiation between maps with and without strings. 

The first step in our implementation of the Canny al- 
gorithm is the filtering of the data. This filtering is in- 
tended to eliminate point source noise. In studies like 
the present one in which we are dealing with simulated 
data which has no noise, the maps can be used without 
smoothing. In order to eliminate point sources, smooth- 
ing should be used in the case of real data. We have 
explored the effects which smoothing has. 

The routine first produces an un-normalized filter map 
F{i,j) (i and j are integer labels running from 1 to Nf 
indexing the pixels of the filter map) 



(16) 



in whichx = i-{Nf-l)/2, y = j-{Nf-l)/2, a = 0.5iV/, 
and Nf is the number of pixels along one direction of the 
filter map, an input parameter in the smoothing routine. 
The normalization is the sum of all of the unnormalized 
weights 



(17) 



Now we can determine the normalized filter map as 



(18) 



Finally, we convolve the input temperature map with the 
filter map to obtain the filtered data map. If M{i,j) 
denotes the initial map, then the filtered map FM{i,j) 
is given by 

FM{i, j) = E - ^' - ' (19) 

k,l=l 

where x and y are determined from k and I as indicated 
below (|16p. Boundary points require special treatment. 
We repeat the points at the boundary. 

The second part of the algorithm involves construct- 
ing the gradient maps. Two arrays are created, the first 
containing the magnitude of the gradient at each pixel 
of the map, the second containing the information about 
the direction of the gradient. 

We compute the lattice derivatives by shifting the tem- 
perature map in each of the eight directions and looking 
for the maximal difference. Instead of using the usual 



lattice gradient which is defined by 

G{x,y) = G{x + e,y)^G{x~e,y)^ 



2e 

G{x,y + e) - G{x,y - e) 
2e 



y, (20) 



(where the grid spacing is e), we shall use the following 
as our definition of the magnitude of the gradient: 



G{x,y) = max|AG| 



(21) 



and the direction is taken to be that for which the max- 
imum value of |AG| occurs. AG is given by the discrete 
gradient 



AG(x,2/) = 



G{n{x,y)) - G{x,y) 
d{x,y) 



(22) 



where n{x, y) denotes the neighboring point to (x, y) and 
d{x,y) is the distance between n{x,y) and {x,y). For 
neighboring points not on the diagonal, the distance be- 
tween the neighboring point and the point {x,y) is 1, for 
points along the diagonal the distance is ^/2. 

The gradient of the image is then stored as an array 
with a dimension containing the two dimensional pixel 
array with the edge strengths and a dimension containing 
the two dimensional pixel array with the edge directions 
as the values. 

The next step is to find the pixels of the map which cor- 
respond to local maxima in the direction of the gradient. 
We first shift the map one unit forwards and backwards 
in each of the four directions (along the coordinate axes 
and along the two diagonals, respectively). For each of 
these four shiftings we focus on points for which the ab- 
solute value of the gradient is greater than at the two 
neighboring points on either side. For all points selected, 
we check if the gradient is in the direction of the shifting. 
Pixel points which pass this test are kept. The others 
are assigned a value of 0. Boundary pixels are stripped 
from the list of local maxima since their magnitudes are 
influenced by boundary effects. In Appendix A we dis- 
cuss the local maxima doubling problem which the code 
has to deal with. 

As mentioned at the beginning of this section, our im- 
plementation of the Canny algorithm makes use of three 
thresholds U and tc. Their values are chosen to corre- 
spond to the amplitude of the KS effect of a single string 
since it is this effect which we want to identify. 

Thus, to set the thresholds for a given value of G/z, we 
first run the code on a number of pure string simulations 
for that value of G/i. We determine the maximal gradi- 
ent (amplitude) Gmax,i for each simulation, and define 
the "string maximal gradient" G„i as the average of the 
individual maxima Gmax,i- 

Returning to our list of candidate local maxima, then 
if the amplitude is larger than tcGm the point is dis- 
carded since it is not due to a typical string gradient. 
For amplitudes in the range 



(23) 
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then the pixel coordinates and gradient directions are 
stored in arrays. Later on, the edge pixel map will be 
constructed by marking all pixel points in the above- 
mentioned coordinate array as 1. These are pixels for 
which the gradient is in the right range to be due to the 
KS effect from a cosmic string. Pixels with a slightly 
lower gradient amplitude might still be due to the KS 
effect from a string with amplitude slightly reduced by 
noise. It is advantageous for the success of the Canny 
algorithm to take into account this possibility. Thus, if 
the amplitude of the local maxima is in the range 

tlGm < A < taGm (24) 

then the coordinates and gradient directions of such pix- 
els are stored in arrays corresponding to pixel points 
marked 1/2. Arrays are created to store the indices of 
the points marked as 1 and 1/2 and separate arrays are 
created to store the directions of the corresponding in- 
dices. 

The next step is the "edgefinder" routine which decides 
whether grid points marked as 1/2 belong to an edge or 
not. Roughly speaking, points marked as 1/2 are consid- 
ered as belonging to an edge if they are direct or indirect 
neighbors to points labelled by 1 and the gradients are 
in a similar direction. To specify the routine, there are 
three choices which can be made for this algorithm, as 
will emerge below. We start from a point labelled 1 (an 
option in the program is to start with a point labelled 
1/2 - the first of the three choices mentioned above)^. 
The program then searches in all eight directions for a 
contiguous point labeled as 1/2 or 1. For each direction, 
we check if it lies in the six directions (or 2 directions) 
perpendicular or near perpendicular (or just perpendic- 
ular - the second of the choices mentioned above) to the 
original point's gradient and check if the gradient is in 
the appropriate parallel or near parallel (or just parallel 
- the third choice) direction to that at the original 1 (or 
1/2). If a point labeled as 1/2 is found, and the direction 
is in one of the allowed six (or two), then we mark the 
coordinates of that point and repeat the search starting 
from that point. We stop if a 1 is found, in which case 
we convert all of the l/2s found on the way to 1. If we 
do not find another point labelled 1/2, we stop and do 
not change the labelling of the points encountered on the 
way. Then, we move on to the next point labelled by 
1 which has not already been covered by the search. In 
the variant of the code in which the search starts with 
points labelled by 1/2, we search for points labelled by 
1/2 or 1 as described above, mark the coordinates of the 
points found, and continue until either a point marked 1 
is found (in which case all of the 1/2 found are converted 
to 1), or else no new point is found, in which case the 



Note that the algorithms starting from 1 or 1/2 are in fact 
different. This is also noted in footnote [7] 



1/2 are considered not to belong to a common edge and 
are not relabelled. 

The output of the "edge- finder" routine is a list of coor- 
dinates of the pixels marked by 1 and of the correspond- 
ing gradient directions. The result of the algorithm can 
be represented by two maps. The first is a map of pixel 
points originally labelled as 1 and 1/2 (different shades), 
the second is an initial "edge pixel map", the map of 
the set of pixels labelled at the end of the "edgefinder" 
routine by a 1. 

The final part of the Canny algorithm is an "edge- 
counting" routine (this part of our method goes beyond 
the standard Canny algorithm). This routine creates a 
second "edge pixel map" and a histogram of number of 
edges as a function of the edge length. Starting with 
the first entry in the list of pixels marked by 1, the rou- 
tine searches in directions perpendicular or near perpen- 
dicular (or exactly perpendicular only) to the gradient 
direction to find other pixels labelled as 1. If such a 
pixel is found, then the routine checks if the gradient di- 
rection of the new pixel is near parallel or parallel (or 
exactly parallel only) to the direction of the gradient at 
the initial pixel. If this is the case, the routine considers 
the new pixel as part of the edge it is following. Once 
no further pixel is found, the routine considers the edge 
to have ended and saves the result. The program then 
moves to the next pixel labelled by 1. Pixels are only 
allowed to be counted in one edge. After all edges are 
found, the lengths are tallied and can be made into a 
histogram. To take into account the fact that a pixel 
may be missing from a string edge due to a large influ- 
ence of the Gaussian component, there is an option for 
the "edge-counting" routine to allow for a gap in an edge 
of a certain length. This skipping length parameter is a 
further input parameter which the user of the algorithm 
has to set. As should also be obvious from the discussion 
in this paragraph, the same choices of "perpendicular ex- 
actly or nearly perpendicular" or "perpendicular exactly" 
and "parallel exactly or nearly parallel" or "exactly par- 
allel" as in the "edge-finding" routine are open to the 
user in this part of the program. 

Let us conclude this section by reminding the reader 
of the various parameters which have to be chosen and 
choices which have to be made in our Canny algorithm 
implementation. First of all, in our "edgefinder" routine 
there are the three thresholds t^^ i« and ti which are 
used in the labelling of the local maxima. Next, in our 
routine to turn 1/2 pixels into 1 pixels, there are three 
choices to be made, first whether one starts with 1 or 
1/2 pixels, second whether one searches for neighboring 
Is and l/2s in direction perpendicular or in directions 
nearly perpendicular or perpendicular to the direction of 
the gradient, and third whether one demands that at the 
selected neighboring site the gradient is parallel or if it is 
nearly parallel or parallel to the direction of the gradient 
at the initial point. In the edge counting algorithm there 
are the two choices analogous to the two choices men- 
tioned at the end of the previous discussion. The user of 
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FIG. 9: Edge maps for the gaussian map shown before. The 
first edge map is for a Gm and thresholds used in the analysis 
of an iV = 10 gaussian map with strings for G/i = 3.5 x 10~^. 
The thresholds are tc — 2, tu = 0.25, ti — 0.1. The second 
edge map is made for the same gaussian map used in the 
analysis of an A'^ = 10 gaussian map with strings for Gfi = 
2 X 10"*. The thresholds are tc = 4, t„ = 0.25, ti = 0.1. 




FIG. 10: The test string edge map. 




FIG. 11; N=l pure string edge map. 



IV. ANALYSIS 



the program has similar choices of "perpendicular exactly 
or nearly perpendicular" versus "perpendicular exactly" 
to the gradient for the directions to search and "paral- 
lel exactly or nearly parallel" versus "parallel exactly" to 
make. Finally, there is the number of points which the 
"edge-counting" routine can skip. 



We tested the capability of our implementation of 
the Canny algorithm to distinguish between tempera- 
ture anisotropy maps which arise from a pure Gaussian 
ACDM model and a model in which there is a contribu- 
tion to the power spectrum coming from a scaling dis- 
tribution of cosmic strings (and with the power of the 
Gaussian noise reduced to maintain the optimal fit to the 
CMB angular anisotropy power spectrum as described in 
Section ini). 
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FIG. 13: = 
noise edge map. 
=0.1. 



2 X 10 unsmoothed string plus gaussian 
N=10. The thresholds are U = 4, f„ = 0.25, 



Since we have applications to small-scale CMB 
anisotropy experiments in mind, we choose a simulation 
box of edge length 10° and with angular resolution of 
1.5 . Both in the ACDM model and in the cosmic string 
model, the actual universe is a realization of a random 
process. Thus, in order to determine whether the Canny 
algorithm can distinguish between the predictions of the 
two theories, we have to run many realizations of the 
models. We ran 50 simulations for both theories. This 
number was chosen since it is the number for which the 



resulting probabilities converge. 

For small values of G/i, we generated each cosmic string 
map independently. For larger values of G/i we generated 
the temperature maps by rescaling the temperature map 
for Gfj, = 6 X 10~®. The value of Gm was scaled accord- 
ingly. The value of the scaling factor a was computed 
for each of the simulations for fixed value of inde- 
pendently. Because of large fluctuations in the various 
random variables (e.g. the number of strings), there re- 
sults a spread for the values of a, and some of the values 
of were negative for large values of G/x. Thus, for large 
values of G/i we used a common value for a^, namely the 
mean of the values of for G/t = 6 x 10~^ and scaled 
according to (see [isj): 



= 1 - 



= 1 - 



iCf)o / Gm 
(GP)o \Gfio 



1 + (ao - 1) 



Gfj_ 
G/io 



25) 



(26) 



where the subscript indicates the reference value of G/i. 

Both the models and the Canny algorithm contain a 
number of parameters. As described in Section |lTl we 
chose the best fit parameters of a ACDM model. For the 
cosmic string distribution, we considered G/i to be the 
free parameter of interest. We fixed the number of strings 
per Hubble volume to be iV = 1 in some simulations and 
N ^ 10 in other simulations, the string segment length 
to be 7_ff^^ with 7 = 1, and the velocity parameter v = 
0.15. 

In this work we did a rough optimization of the var- 
ious parameters which have to be chosen in the Canny 
algorithm. For a fixed value of G/t for which the effect 
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of the strings is statistically significant we varied the pa- 
rameters to find the parameter values which gave the best 
discriminatory power (see below). On this basis, we chose 
the various thresholds. Concerning the optimization of 
the "flags" in the routine, it proved advantageous to start 
with the pixels marked 1 in the edge finding algorithm, 
and to look in directions both perpendicular and near to 
perpendicular, and allowing the gradient at the neigh- 
boring point to be parallel or near to parallel in both the 
edge finding and the edge counting routines 

We ran simulations both with and without smoothing 
of the maps. In the runs with smoothing, the smooth- 
ing was done in the final maps (after adding the string 
maps to the Gaussian maps). Finally, we ran simula- 
tions allowing for skipping of points in the edge counting 
algorithm. 

Discriminating power was quantified using the Fisher 
combined probability test. This test is applied to the two 
output histograms of edge lengths, one from the pure 
Gaussian simulations, the other from the strings plus 
Gaussian maps. For each length I, we are given the 
mean number of edges of that length and the correspond- 
ing standard deviation. Given the two means and corre- 
sponding standard deviations we can apply the t-test to 
compute the probability pi that the two means come from 
the same distribution. The Fisher combined probability 
method then computes a &s follows 



Gmu4d-8. 



k 

1=1 



(27) 



where k is the number of edge lengths being considered, 
and computes the corresponding probability value from 
a distribution with 2fc values. We chose k to be the 
last edgelength to have a nonzero standard deviation in 
either the string or gaussian distributions, whichever is 
smaller. 

The Canny routine makes use of three thresholds. To 
find the cutoff or top threshold we ran a script to find 
the maximum average gradient of the fifty string maps 
and then compared the highest gradient to the highest 
average gradient. This threshold needs to increase as the 



We ran the Canny algorithm on two different test maps to de- 
velop an intuition on how to optimize the flags. For both test 
maps, more edges were found when the number of allowed direc- 
tions and gradients were greater. In one test map more edges 
were found when the algorithm started with points marked as 
1/2 and in the other algorithm more edges were found when the 
algorithm started with points marked as 1 for some values of the 
flags. Since the algorithm is much slower when it starts with 
points marked as 1/2, in this paper we always start with points 
marked as 1. Longer edges were found when skipping points was 
allowed. This is consistent with the findings we have for the 
averages of the full 50 maps of gaussian temperatures and the 
combined maps with strings. Since skipping reduces the number 
of short edges, where the standard deviation is smaller compared 
to the mean, and thus increases the p- value, it increases the prob- 
ability that the maps look the same. 
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FIG. 15: Histogram of the edge length distribution of the un- 
smoothed pure Gaussian maps and the strings plus smoothed 
Gaussian maps for a value of G^. = 4 x 10~*. The number of 
edgelengths for strings plus gaussian maps are denoted by tri- 
angles and the squares denote the number of edgelengths for 
the pure gaussian maps. We plot the means and the standard 
deviation as the error. Note the systematically larger number 
of edges for maps with strings for short edgelengths. 



Gp, decreases because otherwise one throws out virtually 
all of the edges, since as G/z decreases Gm decreases. One 
needs to keep some edges with values greater than the 
maximum gradient in any string map because the signal 
from the gaussian map may add in the same location to 
the signal from the string map. To find the upper and 
lower thresholds we ran the Canny algorithm for various 
thresholds on a pure string map and determined by eye 
if enough strings could be seen but not so many that the 
overlap of the edges breaks apart the long string edges. 
In particular, for = 10 if the thresholds are too small 
the images are saturated with short edges which break up 
the long edges. Thus, the thresholds need to be higher for 

= 10 than for iV = 1. We leave a detailed optimization 
of the thresholds to future work. 

As a test of the code, we ran the analysis algorithm on 
two independent sets of fifty Gaussian CMB maps and 
computed the probability that the two sets come from 
the same Gaussian ensemble. The resulting probability 
was typically of the order of 0.7, i.e. within the la error 
range. 

Our simulations show that maps originating from cos- 
mic strings have a larger number of edges than the corre- 
sponding pure Gaussian maps. This is a consequence of 
the presence of the Kaiser-Stebbins edges, the signal we 
are looking for. Figure [T51 shows a comparison of the two 
histograms for the value G/i ~ 4 x 10^^. Since the Gaus- 
sian noise cuts up the long edges, the difference in edge 
lengths is statistically significant only for short edges. 
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Table 1 summarizes our results obtained by apply- 
ing the Fisher combined probability method to the his- 
tograms of string edge lengths. The results are for 
N = 10 (10 strings per Hubble volume per Hubble 
time), for un-smoothed maps, without skipping in the 
edge drawing and edge counting routines. The statisti- 
cal analysis includes edgelength 1 data. The first column 
gives the value of G/i, the second the value of the top 
threshold (the other two thresholds are held fixed at the 
values tu = 0.25 and ti =0.1), the third gives the prob- 
ability that the histograms of the string and Gaussian 
maps come from the same Gaussian distribution. The 
fourth value gives the maximal edge length used in the 
analysis, and the last column lists the value of a. We 
see that the cosmic string signal can be detected ^ at a 
three sigma level down to a value — 2 x 10~^ which 
is almost two orders of magnitude better than the cur- 
rent bounds obtained by direct searches, and one order of 
magnitude better than the limits on G/i coming from the 
CMB angular power spectrum constraints. These results 
are for unsmoothed maps. 

As can be expected from looking at the histograms of 
[T5l the sensitivity of the algorithm decreases slightly if 
edges of length 1 are excluded from the analysis. The 
results excluding edges of length 1 are shown in Table 2. 

Since smoothing dramatically reduces the number of 
short edges, and since the strength of the signal of our 
analysis comes from the number of short edges, it turns 
out that smoothing significantly weakens the ability of 
the algorithm to pick out strings. The results of our 

* More conservatively, we should say that the difference between 
the maps with cosmic strings and pure Gaussian maps can be 



analysis applied to smoothed maps (smoothing length 3) 
are given in Table 3. The loss in discriminatory power 
is about a factor of 5. This is a serious concern when 
considering applications of our algorithm to real data. 

In Table 4 we present the results for iV = 1. Since 
there are an order of magnitude less string edges in this 
case, the limit on Gfi which can be obtained is slightly 
weaker (about a factor of 3 weaker). To put this result 
into context, it is important to point out that limits on 
Gfi from other studies implicitly or explicitly use N — 
10. The limits on G/i from matching the angular CMB 
power spectrum would be weaker by about one order of 
magnitude. 

The results of the previous tables were obtained from 
an algorithm which did not have any skipping in the edge 
counting routine. We allowed skipping of two points. 
Skipping of two points leads to a larger number of long 
edges. Skipping by four points leads to edge broadening, 
an unwanted feature. Hence, we only considered skip- 
ping two points. Including skipping turns out to reduce 
the number of short edges more significantly than it in- 
creases the number of longer edges. Thus, the power of 
our routine to discriminate between maps with and with- 
out strings slightly decreases when introducing skipping 
(which is the opposite of what we initially expected) . The 
results are indicated in the following table. However, it is 
possible that further optimization of the routine parame- 
ters would reverse the results concerning the effectiveness 
of skipping. 

detected. 



Table 1: Unsmoothed String Map versus Gaussian Map N=10 
Minimum Edge Length: 1 
tu = 0.25, ti = 0.1, num. skipped points=0 



Gmu 




probability 


max. edgelength 


mean a 


l.OOOOOOOe-08 


7.00000 


0.44757170 


11.000000 


0.99967664 


2.0000000e-08 


4.00000 


0.0020638683 


14.000000 


0.99862559 


4.0000000e-08 


3.00000 


0.0000000 


21.000000 


0.99447207 


6.0000000e-08 


3.00000 





33 


0.98754092 



Table 2: Unsmoothed String Map versus Gaussian Map: N=10 
Minimum Edge Length: 2 
tu = 0.25, ti = 0.1, skipped points=0 



Gmu 
l.OOOOOOOe-08 


7.00000 


probability 
0.56022527 


max. edgelength 
11.000000 


mean a 
0.99967664 


2.0000000e-08 


4.00000 


0.087458606 


14.000000 


0.99862559 


4.0000000e-08 


3.00000 


4.8140603e-ll 


21.000000 


0.99447207 


6.0000000e-08 


3.00000 


0.0000000 


33.000000 


0.98754092 



V. CONCLUSIONS AND DISCUSSION 

We have developed a new program to search for cos- 
mic strings in CMB anisotropy maps making use of the 



Canny algorithm, and have tested it with simulated data 
corresponding to maps with specifications corresponding 
to those of current ground-based CMB experiments. The 
code contains a number of optimization parameters, and 
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Table 3 
Minimum Edge Length: 1 
Smoothed String Map versus Gaussian Map N=10 
tu = 0.25, ti = 0.1, num. skipped points=0 



Gmu 


tc 


probabihty 


max. edgelength 


mean a 


6.0000000e-08 


3.00000 


0.10701445 


42.000000 


0.98754092 


8.0000000e-08 


3.00000 


0.059179809 


46.000000 


0.97932855 


9.0000000e-08 


3.00000 


0.00012533417 


46.000000 


0.97174247 


l.OOOOOOOe-07 


3.00000 


2.4759972e-ll 


52.000000 


0.96499435 


1.5000000e-07 


2.50000 


0.0000000 


59.000000 


0.91936485 



Table 4 
Minimum Edge Length: 1 
Unsmoothed String Map versus Gaussian Map: N=l 
tu ~ 0.03, ti = 0.005, skipped points=0 



Gmu 




probability 


max. edgelength 


mean a 


2.0000000e-08 


4.00000 


0.43985131 


7.0000000 


0.99986506 


4.0000000e-08 


3.00000 


0.76437603 


11.000000 


0.99937157 


6.0000000e-08 


3.00000 


0.0019923380 


19.000000 


0.99901436 


8.0000000e-08 


3.00000 


8.7787192e-09 


28.000000 


0.99816180 


9.0000000e-08 


3.00000 


3.8369308e-13 


32.000000 


0.99757307 



we have discussed the role of these parameters. 

Our results (based on a rough optimization of the pa- 
rameters) show that our algorithm has the potential to 
improve the bounds on the cosmic string tension from di- 
rect CMB observations by up to two orders of magnitude 
compared to existing limits, and by one order of magni- 
tude from other more indirect limits. The limiting value 
of decreases as N (the number of strings per Hubble 
volume per Hubble time) increases. For cosmic super- 
strings the intercommutation probability may be much 
lower than for Abelian field theory strings, thus leading 
to a larger value of N [HB]- Hence, our method may be 
able to set more stringent limits of the value of G/i in the 
case of cosmic superstrings. ^ 

A possible concern with our analysis is that our toy 
model of representing the effects of a string on the mi- 
crowave sky introduces artificial edges. Although the 
presence of the central edge in the CMB pattern of a 
string is clearly physical, the edges perpendicular to this 
central edge are due to the toy model which considers 
finite length string segments rather than a network of 
infinite strings. Although the temperature perturbation 
clearly goes to zero at a perpendicular distance t from 
a string pres ent at time t, the sharpness of the change 
given by [2g| may be modulated by interactions between 
the time of formation of the string network and the time 
t Thus, one of our simulations with N strings per 
Hubble volume per Hubble time may contain as many 
edges as a "real" string map with 4A^ strings per Hubble 



While this manuscript was being prepared for submission, a 
preprint [59ll appeared in which a novel method different from 
ours was proposed to search for cosmic strings. With this new 
method, improved Umits on the cosmic string tension also appear 
to be possible. 

We thank the Referee for raising this issue. 



volume per Hubble time. As can be seen by compar- 
ing Tables 4 and 5, the dependence of the limits on G/i 
on N is not large. The limit on G/x appears to scale as 
N~^. Thus, the limit on obtained from our analysis 
might be a factor of 3 too optimistic due to the larger 
number of edges. This point could also be addressed by 
a Canny algorithm analysis of more sophisticated string 
simulations, which we shall pursue as a follow-up study. 

The results presented here are based on comparing sim- 
ulations with and without cosmic strings of a known value 
of Gfi. However, we do not know the value of Gfj,, and 
the goal of observations is to determine the value or to 
set limits on it. Hence, the way we plan to apply our 
algorithm to real data is to compare the edge histograms 
based on real data with those of a cosmic string plus 
Gaussian map with fixed G/i, and varying Gfi, starting 
from the current limit for this quantity. Limits on Gfi 
can be obtained if the histogram derived from the data is 
statistically different from that of a Gaussian plus strings 
map for the corresponding value of the normalized string 
tension. A candidate detection of cosmic strings would 
require the histogram derived from the real data to be 
statistically indistinguishable from that of a Gaussian 
plus strings map, and at the same time different from 
that of a pure Gaussian map. If there are indeed cos- 
mic strings in the sky with a value of G/i for which our 
analysis gives a probability p that the histograms with 
and without strings are from the same distribution, then 
the probability that our statistical analysis would give a 
positive detection of these strings will be 1 — p. 

When applying the algorithm to real data, a serious 
concern is systematic errors introduced by the observing 
strategy - specifically lines in the maps due to the scan- 
ning strategy. In specific experiments such as the South 
Pole Telescope (SPT) experiment, lines introduced by 
the scanning will be in one specific direction, and thus 
the algorithm can be rendered immune to this effect by 
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Table 5 
Minimum Edge Length: 1 
Unsmoothed String Map versus Gaussian Map:N=10 
tu ~ 0.25, ti — 0.1, skipped points=2 



Gmu 




probability 


max. edgelength 


mean a 


l.OOOOOOOe-08 


7.00000 


0.43783402 


41.000000 


0.99967664 


2.0000000e-08 


4.00000 


0.30944903 


55.000000 


0.99862559 


4.0000000e-08 


3.00000 


5.4179907e-08 


139.00000 


0.99447207 



considering only edges which are not parallel to the di- 
rection of the scanning stripes. Foreground and instru- 
mental noise are other important problems. Smoothing 
of the maps is supposed to reduce the latter problem, and 
the study of |48j confirms that instrumental noise such 
as that anticipated in the SPT telescope will not have a 
big effect on the power of the algorithm. 

The next step of our research program will be to apply 
our code to existing data. The main challenge will be to 
take care of the systematic errors contained in the data. 
We plan to compare data for the actual sky map to Gaus- 
sian and Gaussian plus strings simulations. In order to 
put the data in a format for which a statistical analysis 
is possible, we will divide the sky maps into sub-maps 
of equal size, and treat each sub-map as a different real- 
ization of the sky. In other words, we make the ergodic 
hypothesis and identify ensemble and spatial averaging. 
Thus, if we want to compare squares of edge length 10*^ in 
the sky, then we need a total survey area 50 times larger 
in order to be able to have 50 independent data sets for 
the sky. We also plan to apply our code to more realistic 
cosmic string simulations, such as those of [53 |. 
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VI. APPENDIX A 

In this Appendix we describe the maxima doubling 
problem which arises in our implementation of the Canny 
algorithm. Let us illustrate the problem in terms of an 
example. 

If the temperature grid points are 
6 8 10 20 

3 5 12 14 (28) 

4 6 9 11 

and the points are labeled starting from from left to 
right top to bottom, then the point at grid point 5 (tem- 
perature 5) is a local maximum with gradient 7 along the 
direction. The grid point 6 (temperature 12) is a local 
maximum along the 1 direction with gradient 8. There 
aren't any doubles in this case. 

However, if we had 

6 8 10 15 

3 5 12 14 (29) 

4 6 9 11 

then the grid point 5 (temperature 5) is a local maximum 
along the direction with a gradient of 7 and the grid 
point 6 (temperature 12) is a local maximum along the 
4 direction with gradient 7. We use a sorting routine to 
remove the first occurrence of a gradient of 7 (index 5) 
as we want there to be only one local maximum at this 
point. 



[1] A. C. Davis and T. W. B. Kibble, "Fundamental 
cosmic strings," Con temp. Phys. 46, 313 (2005) 
arXiv: hep-th/0505050'; 

M. Sakellariadou, "Cosmic Superstrings," 

arXi^:0802.3379 [hep-th]. 
[2] R. Jeannerot, "A Supersymmetric SO (10) Model with 
Inflation and Cosmic Strings," Phys. Rev. D 53, 5426 
(1996) ^arXiv:h ep-ph/9509365 ; 

R. Jeannerot, J. Rocher and M. Sakellariadou, "How 
generic is cosmic string formation in SUSY GUTs," Phys. 
Rev. D 68, 103514 (2003) arXiv:hep-ph/0308134 . 
[3] S. Sarangi and S. H. H. Tye, "Cosmic string production 



towards the end of brane inflation," Phys. Lett. B 536, 
185 (2002) [arXiv:hep -th/0204074 . 
[4] E. Witten, "Cosmic Superstrings," Phys. Lett. B 153, 
243 (1985). 

[5] E. J. Copeland, R. C. Myers and J. Polchinski, 
"Cosmic F- and D-strings," JHEP 0406, 013 (2004) 
arXiv;hep-th/0312067 . 

[6] J. Khoury, B. A. Ovrut, P. J. Steinhardt and N. Turok, 
"The ekpyrotic universe: Colliding branes and the origin 
of the hot big bang," Phys. Rev. D 64, 123522 (2001) 
fM-Xiv:hep-th/0103239 . 

[7] R. H. Brandenberger and C. Vafa, "Superstrings In The 



18 



Early Universe," Nucl. Phys. B 316, 391 (1989). 
[8] A. Nayeri, R. H. Brandenberger and C. Vafa, "Producing 
a scale-invariant spectrum of perturbations in a Hage- 
dorn phase of string cosmology," Phys. Rev. Lett. 97, 
021302 (2006) arXiv:hep-th/0511140 ; 
R. H. Brandenberger, A. Nayeri, S. P. Patil and 
C. Vafa, "String gas cosmology and structure formation," 
'arXiv:hep-th/0608121 
[9] T. W. B. Kibble, "Phase Transitions In The Early Uni- 
verse," Acta Phys. Polon. B 13, 723 (1982); 
T. W. B. Kibble, "Some Implications Of A Cosmological 
Phase Transition," Phys. Rept. 67, 183 (1980). 

[10] A. Vilenkin and E.P.S. Shellard; Cosmic Strings and 
Other Topological Defects, (Cambridge Univ. Press, 
Cambridge, 1994). 

[11] M. B. Hindmarsh and T. W. Kibble, "Cos- 
mic strings," Rept. Prog. Phys. 58, 477 (1995) 
,arXiv :hep-ph/9411342^. 

[12] R. H. Brandenberger, "Topological defects and struc- 
ture formation," Int. J. Mod. Phys. A 9, 2117 (1994) 
arXiv:astro-ph/9310041 . 

[13] T. W. B. Kibble, "Topology Of Cosmic Domains And 
Strings," J. Phys. A 9, 1387 (1976); 

Y. B. Zeldovich, "Cosmological fluctuations produced 
near a singularity," Mon. Not. Roy. Astron. Soc. 192, 
663 (1980); 

A. Vilenkin, "Cosmological Density Fluctuations Pro- 
duced By Vacuum Strings," Phys. Rev. Lett. 46, 1169 
(1981) [Erratum-ibid. 46, 1496 (1981)]. 

[14] A. Vilenkin, "Cosmic Strings And Domain Walls," Phys. 
Rept. 121, 263 (1985). 

[15] A. Albrecht and N. Turok, "Evolution Of Cosmic 
Strings," Phys. Rev. Lett. 54, 1868 (1985). 

[16] D. P. Bennett and F. R. Bouchet, "Evidence For A Scal- 
ing Solution In Cosmic String Evolution," Phys. Rev. 
Lett. 60, 257 (1988). 

[17] B. Allen and E. P. S. Shellard, "Cosmic String Evolu- 
tion: A Numerical Simulation," Phys. Rev. Lett. 64, 119 
(1990). 

[18] G. Vincent, N. D. Antunes and M. Hindmarsh, "Nu- 
merical simulations of string networks in the Abelian- 
Higgs model," Phys. Rev. Lett. 80, 2277 (1998) 
arXiv:hep-ph/9708427 . 

[19] V. Vanchurin, K. D. Olum and A. Vilenkin, "Scaling of 
cosmic string loops," Phys. Rev. D 74, 063527 (2006) 
arXiv:gr-qc/0511159 . 

[20] C. J. A. Martins and E. P. S. Shellard, "Frac- 
tal properties and small-scale structure of cosmic 
string network s," Phys . Rev. D 73, 043515 (2006) 
arXiv: ast ro-ph/05 1 1 7921 . 

[21] C. Ringeval, M. Sakellariadou and F. Bouchet, "Cosmo- 
logical evolution of cosmic string loops," JCAP 0702, 
023 (2007) arXiv:astro-ph/0511646 . 

[22] J. N. Moore and E. P. S. Shellard , "On the_evolution_o f 
abelian-Higgs string networks," |arXiv:hep-ph/9808336 

[23] J. Polchinski, "Cosmic String Loops and Gravitational 
Radiation," arXiv:0707.0888 [astro-ph]. 

[24] H. B. Nielsen and P. Olesen, "Vortex-Line Models For 
Dual Strings," Nucl. Phys. B 61, 45 (1973). 

[25] N. Turok and R. H. Brandenberger, "Cosmic Strings And 
The Formation Of Galaxies And Clusters Of Galaxies," 
Phys. Rev. D 33, 2175 (1986); 

H. Sato, "Galaxy Formation by Cosmic Strings," Prog. 
Theor. Phys. 75, 1342 (1986); 



A. Stebbins, "Cosmic Strings and Cold Matter", Ap. J. 

(Lett.) 303, L21 (1986). 
[26] N. Kaiser and A. Stebbins, "Microwave Anisotropy Due 

To Cosmic Strings," Nature 310, 391 (1984). 
[27] A. Vilenkin, "Gravitational Field Of Vacuum Domain 

Walls And Strings," Phys. Rev. D 23, 852 (1981). 
[28] J. C. R. Magueijo, "Inborn metric of cosmic strings," 

Phys. Rev. D 46, 1368 (1992). 
[29] A. Kosowsky [the ACT Collaboration], "The Atacama 

Cosmology Telescope Project: A Progress Report," New 

Astron. Rev. 50, 969 (2006) arXiv:astro-ph/0608549 . 
[30] J. E. Ruhl et al. [The SPT Collaboration], "The South 

Pole Telescope," Proc. SPIE Int. Soc. Opt. Eng. 5498, 

11 (2004) arXiv:astro-ph/0411122 . 
[31] F. Villa et al. [The Planck Collaboration], "The 

Planck Telescope," AIP Conf. Proc. 616, 224 (2002) 
arXiv:a stro-ph/0112173 ; 

J. Tauber (on behalf of ESA and the Planck Science Col- 
laboration), "The Planck Mission", Advances in Space 
Research (2004). 

[32] R. K. Sachs and A. M. Wolfe, "Perturbations of a cos- 
mological model and angular variations of the microwave 
background," Astrophys. J. 147, 73 (1967). 

[33] A. Albrecht, D. Coulson, P. Ferreira and J. Magueijo, 
"Causality and the microwave background," Phys. Rev. 
Lett. 76, 1413 (1996) (arXiv:astro-ph/9505030 . 

[34] L. Perivolaropoulos, "Spectral Analysis Of Microwave 
Background Perturbations Induced By Cosmic Strings," 
Astrophys. J. 451, 429 (1995) arXiv:astro-ph/9402024; . 

[35] J. Magueijo, A. Albrecht, D. Coulson and P. Ferreira, 
"Doppler peaks from active perturbations," Phys. Rev. 
Lett. 76, 2617 (1996) arXiv:astro-ph/9511042 . 

[36] U. L. Pen, U. Seljak and N. Turok, "Power 
spectra in global defect theories of cosmic struc- 
ture formation, " Phys. Rev. Lett. 79, 1611 (1997) 
[arXiv:astro-ph/9 704 165 . 

[37] R. A. Sunyaev and Y. B. Zeldovich, "Small scale fluc- 
tuations of relic radiation," Astrophys. Space Sci. 7, 3 
(1970). 

[38] P. D. Mauskopf et al. [Boomerang Collaboration], "Mea- 
surement of a Peak in the Cosmic Microwave Back- 
ground Power Spectrum from the North American test 
flight of BOOMERANG," Astrophys. J. 536, L59 (2000) 
arXiv:astro-ph/9911444 . 

[39] C. L. Bennett et al., "First Year Wilkinson Microwave 
Anisotropy Probe (WMAP) Observations: Preliminary 
Maps and Basic Results," Astrophys. J. Suppl. 148, 1 
(2003) arXiv:astro-ph/0302207 . 

[40] L. Pogosian, S. H. H. Tye, I. Wasserman and M. Wyman, 
"Observational constraints on cosmic string pro- 
duction during brane inflation," Phys. Rev. D 68, 
023506 (2003) [Erratum-ibid. D 73, 089904 (2006)] 
arXiv:hep-th/0304188 ; 

M. Wyman, L. Pogosian and I. Wasserman, "Bounds on 
cosmic strings from WMAP and SDSS," Phys. Rev. D 
72, 023513 (2005) [Erratum-ibid. D 73, 089905 (2006)] 
arXiv:astro-ph/0503364 ; 

A. A. Fraisse, "Limits on Defects Formation 
and Hybrid Inflationary Models with Three- Year 
WMAP Observations," JCAP 0703, 008 (2007) 
arXiv:astro-ph/0603589 ; 

U. Seljak, A. Slosar and P. McDonald, "Cosmological 
parameters from combining the Lyman-alpha forest with 
CMB, galaxy clustering and SN constraints," JCAP 



19 



[41] 



[42] 



[43] 



[44] 



0610, 014 (2006) 'arXiv:astro-ph/0604335"; 
N. Bevis, M. Hindmarsh, M. Kunz and J. Urrestilla, 
"CMB power spectrum contribution from cosmic 
strings using field-evolution simulations of the Abelian 
Higgs mo del," Phys. Rev. D 75, 065015 (2007) 
larXiv:astro-ph/0605018 ; 

N. Bevis, M. Hindmarsh, M. Kunz and J. Urrestilla, 
"Fitting CMB data with cosmic strings and inflation," 
arXiv:astro-ph/0702223 

R. A. Battye, B. Garbrecht and A. Moss, "Constraints 
on supersymmetric models of hybrid inflation," JCAP 
0609, 007 (2006) [arXiv:astro-ph/0607339| ; 
R. A. Battye, B. Garbrecht, A. Moss and H. Stoica, 
"Constraints on Brane Inflation and Cosmic Strings," 
JCAP 0801, 020 (2008) arXiv:0710.1541 [astro-ph]]. 
F. R. Bouchet and D. P. Bennett, "Does The Millisecond 
Pulsar Constrain Cosmic Strings?," Phys. Rev. D 41, 720 
(1990). 

V. M. Kaspi, J. H. Taylor and M. F. Ryba, "High - pre- 
cision timing of millisecond pulsars. 3: Long - term mon- 
itoring of PSRs B1855+09 and B1937+21," Astrophys. 
J. 428, 713 (1994); 

S. E. Thorsett and R. J. Dewey, "Pulsar timing limits 
on very low frequency stochastic gravitational radiation," 
Phys. Rev. D 53, 3468 (1996); 

M. P. McHugh, G. Zalamansky, F. Vernotte and 
E. Lantz, "Pulsar Timing And The Upper Limits On A 
Gravitational Wave Background: A Bayesian Approach," 
Phys. Rev. D 54, 5993 (1996); 

A. N. Lommen, "New limits on gravitational radiation 
using pulsars," arXiv:astro-ph/0208572 
'Upper 



F. A. Jenet et al. 



on the low- 



bounds 

frequency stochastic gravitational wave background 
from pulsar timing observations: Current limits and 
future prospects," Astrophys. J. 653, 1571 (2006) 
[arXiv:astro-ph/0609013 . 

R. Moessner, L. Perivolaropoulos and R. H. Branden- 
berger, "A Cosmic string specific signature on the cosmic 
microwave background," Astrophys. J. 425, 365 (1994) 
arXiv:astro-ph/9310001 . 

A. S. Lo and E. L. Wright, "Signatures of cos- 
mic strings in the cosmic microwave background," 



arXiv:astro -ph/0503120 
[45] E. Jeong and G. 



[46] 



F. Smoot, "Search for cosmic strings 
in CMB anisotropies," Astrophys. J. 624, 21 (2005) 
,arXiv:astro-ph/0406432 ; 

E. Jeong and G. F. Smoot, "The Vahdity of the Cos- 
mic String Pattern Search with the Cosmic Microwave 
Background," arXiv:astro-ph/0612706 
S. Amsel, J. Berger and R. H. Brandenberger, "Detecting 



Cosmic Strings in the CMB with the Canny Algorithm," 
JCAP 0804, 015 (2008) arXiv:0709.0982 [astro-ph]]. 
[47] J. Canny, "A computational approach to edge detection" , 
IEEE Trans. Pattern Analysis and Machine Intelligence 
8, 679 (1986). 

[48] A. Stewart and R. Brandenberger, "Edge Detec- 
tion, Cosmic Strings and the South Pole Telescope," 
arXET0809.0865 [astro-ph]. 

[49] M. J. White, J. E. Carlstrom and M. Dragovan, "In- 
terferometric Observation of Cosmic Microwave Back- 
ground Anisotropies," Astrophys. J. 514, 12 (1999) 
[arXiv:astro -ph/9712195 . 

[50] A. Lewis, A. Challinor and A. Lasenby, "EfBcient Com- 
putation of CMB anisotropies in closed FRW models," 
Astrophys. J. 538, 473 (2000) arXiv:astro-ph/9911177; . 

[51] U. Seljak and M. Zaldarriaga, "A Line of Sight Approach 
to Cosmic Microwave Background Anisotropies," Astro- 
phys. J. 469, 437 (1996) arXiv:astro-ph/9603033 . 

[52] C. 1. Kuo et al. [ACBAR coUaborationTT^ "High 
Resolution Observations of the CMB Power Spec- 
trum with ACBAR," Astrophys. J. 600, 32 (2004) 
arXiv:astro-ph/0212289 . 

[53] C. L. Reichardt et al, "High resolution CMB power 
spectrum from the complete ACBAR data set," 
,arXiv:0801.1491 [astro-ph]. 

[54] A. A. Fraisse, C. Ringeval, D. N. Spergel and 
F. R. Bouchet, "Small-Angle CMB Tempera- 
ture Anisotropies Induced by Cosmic Strings," 
arXiv:0708.1162 [astro-ph]. 

[55] L. Perivolaropoulos, "COBE versus cosmic strings: An 
Analytical model," Phys. Lett. B 298, 305 (1993) 
jarXiv:hep-ph/9208247 ; 

L. Perivolaropoulos, "Statistics of microwave fluctuations 
induced by topological defects," Phys. Rev. D 48, 1530 
(1993) arXiv:hep-ph/9212228 . 

[56] J. Polchinski, "Introduction to cosmic F- and D-strings," 
arXiv:hep-th/0412244 

[57] J. H. Traschen, N. Turok and R. H. Brandenberger, "Mi- 
crowave Anisotropies from Cosmic Strings," Phys. Rev. 
D 34, 919 (1986). 

[58] E. Komatsu et al. [WMAP Collaboration], "Five- 
Year Wilkinson Microwave Anisotropy Probe Obser- 
vations:Cosmological Interpretation," larXiv:0803.0547l 
[astro-ph]. 

[59] D. K. Hammond, Y. Wiaux and P. Vandergheynst, 
"Wavelet domain Bayesian denoising of string signal 
in the cosmic microwave background," arXiv:0811.1267 
[astro-ph]. 



